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Abstract 

We solve problem x proposed by O. Oloa, AMM xxx 2012 119? (to ap- 
pear), p. yyy for certain definite logarithmic integrals. A number of generating 

functions are developed with certain coefficients p n , and some extensions are 

(n) 

presented. The explicit relation of p n to Norlund numbers B n is discussed. 
Certain inequalities are conjectured for the {p n } sequence of coefficients, in- 
cluding its convexity, and an upper bound is demonstrated. It is shown that 
p n values may be used to express the Stieltjes constants for the Hurwitz and 
Riemann zeta functions, as well as values of these zeta functions at integer 
argument. Other summations with the p n coefficients are presented. 
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Solution of problem xxx 

0. Oloa has proposed the following problem in the Amer. Math. Monthly 119 
(?), yyy (2012?). 
(a) Prove the formula 

1 / 1 1 V , , ,„ s 3 



o 



+ dx = In(27r) - -. (1.1) 



In x 1 — x I 2 



(b) If a > 0, find a closed form expression for 



1 * i i ^ 2 



In x 1 — x / 

N.B. Although the statement (b) is restricted to o > 0, this part may be extended 
to Re a > — 1, and our explicit expressions reflect this fact. 

We first prove (1.1) and answer part (b), and then present several extensions. 

Proof. Let (a) n = T(a + n)/T(a) be the Pochhammer symbol, where T is the 
Gamma function, = T'/T be the digamma function, 7 = — ^(1) be the Euler 
constant, and p F q the generalized hypergeometric function (e.g., [3[1]). We introduce 
the positive constants (e.g., [3], Proposition 11, [3], Proposition 5, [5J, Proposition 2) 

If 1 , x , (-l) n+1 >As(n,fc) , . 

n\ J n\ k + 1 

where s(/c, t) is the Stirling number of the first kind. The first few values of these are 
p 2 = 1/2, p 3 = 1/12, p 4 = 1/24, p 5 = 19/720, and p 6 = 3/160. These constants enter 

the generating function 

00 

T.P^' 1 = ~ + Tn v 1*1 <L ( L3 ) 

*— ' z ln(l — z) 

n=l v ' 
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Multiplying (1.3) by ln(l — z) and manipulating series, one finds the recursion relation 

n-l 



1 



Pn+l = 



1 4^ (n- j 



n + 



j+iy 



n > 1. 



We then have 
Lemma 1. 



1 x 

ln 2 (1 _ s = + X K+3 - ™Pn+2K l + P3 + 

^ ' n=l 

This follows from the derivative expression 



\z\ < 1. 



oo 11 

J](n - l)p n+1 z n - 2 = --J + 2 

^ ^ (1 - *) In (1 - *) 



As a consequence, we obtain 



Then for (a), 



1 1 

+ 



o KM 1 -*) z 



dz 



1 1 

+ 



In z 1 — z 



dz 



l ji + ]D( n + 3 ^+ 3 - n Pn+2\{l ~ zT^j dz 
x 

+ y^[( n + 3 K+3 - np n+2 ] — — . 



1 

4 



;i-4) 



;i.5) 



/ i i\ 2 i 00 

( ln(1 _ x + -) = i + E[( n + 3 W-n Pn+2 ]^. (1.6) 

\ V / / n=l 



(1.7) 



The integral representation of (1.2) may now be inserted, and the sums rewritten in 
terms of binomial coefficients. For instance, we have 

2 



2^n + l(n + 2)\ { X)n+2 

n=l v 7 



^ ] (n + l)\n + 2j ^ ] 

n=l v ' x ' n=l 



1 + 



(n + 1) 



x 

n + 2 



-(2 - x)(x - 1) - x[l - 27 + x - 2ip{x + 1)]. 



(1.8) 



Noting that ip(x + l)c?x = 0, we have 



ln(l — z) z 



- dz=- + 



x 3 
7 — — — 27X + -x 2 — 2xip(x + 1) 



dx 



ln(2vr) 



(1.9) 



In the last step, we integrated by parts (e.g., [2]), 



xi/j(x+l)dx 



[ lnT(x+l)dx = - [ [lnx+\nT{x)}dx = l--ln(27r). (1.10) 
Jo Jo 2 



For (b), we let B(x,y) = T(x)r(y)/T(x + y) be the Beta function. Then 

x a dx 



1 ' 1 1 ^ 2 



In x 1 — x 



~~ J Q |i + XjK n + 3 ^™+ 3 ~ ^ +2 ](1 - z)"J 
77-— pr + Y\[(n + 3K +3 - n Pn+2 ]B(n + l,a + l). (1.11) 



Lemma 2. For Re y > 0, 

^p n+3 -B(ri + l,y) = -- 



n=l 



_L _ l + 2y - ln(2vr) + 21nr(y) + (1 - 2y)^(y) 
6y 



(1.12) 



Proof. From (1.2) we have 



00 1 r 1 

J2p n+3 B{n + l,y) = -— x{x-l)[ 3 F 2 {l,l,2-x;3,y + l;l)-l]dx. (1.13) 



We now use the identity 



(l)i = 2 (l)i _ (2)j 
(3)i (2), (3), : 



(1.14) 



being a special case of 



K- =(a + 1) _J a + 1 ^ 



to obtain 



(a + 2), > + !); (a + 2),' 



E 

3=0 



(1), (2), 



(2), (3), 



(2 - (1), 



2y 

-' -[l _ x - (x + 3/ - 1)V(S/) + (a: + 2/ - + ?/ + 1)]. (1.15) 



x(x — 1) 

Carrying out the integration of (1.13) gives the Lemma. □ 
Then by Proposition 2 in the next section with y = a + 1 we obtain 

1 ■ / 1 1 \* „ ,, _ „ _ 3 



+ - x°dx = (o-+l)ln(f7+l)-2(7+(7V'((7+l)-21nr((7+l)+ln(27r)--. 

Ill CC .L X J Zi 

(1.16) 

Remarks. As is apparent from (1.11) and (1.16), Re a = —1 is the 'critical line' 
for divergence of the integral. 

We have the following 
Corollary 1. For n > 1, 

n+l 

(n + 3)p n+3 - np n+2 = ^p fc+ ip n _ fc+3 . (1.17) 

fc=i 

This follows from multiplication of series, using (1.6), 
/ 1 1\ 2 1 °° 

V ln(l-*) + ~z) = l + J2tt n + ^P n +3- n P n +^ zn 
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oo n+1 

= J2J2 pk + iPn - k + 3zn - ( L18 ) 

n=0 k=l 

Is there a combinatorial interpretation of identity (1.13)? 

The Appendix generalizes (1.14) and the following identity for ratios of Pochham- 
mer symbols. 

Extensions 

We may proceed similarly as above, and find for instance 

1 1 °° 

In 1 — z) I ^— ' 

v y n=l 

giving, along with (1.3) and (1.4), 

1\ 3 1 3 3 1 



X^X-z)^ z) ln 3 (l-^) + zln 2 (l-^) + ^ 2 ln(l-^) + z s 
- Y}(n + 4)(n + 5K +4 - (ra + l)(2ra + 7)p n+3 + n 2 pn+ 2 ]^ 



2 

n=4 



Then 



1 z z 2 133 o . , 

+ 8 + 16 " 24 + 4320' (2 ' 2) 



i 1 \ 3 r 1 / i i x 3 



+ -) dz= - — + cte 



o \ln(l — z) z) J Q \\nz 1 — z t 

\Yyin + 4)(n + 5K +4 - (n + l)(2n + 7)p n+3 + ^WJt^t + ^| 

n=4 ^ ' 

= -^ + 61nA (2.3) 



6 



wherein A is Glaisher's constant, such that hi A = — [C(— 1) + C'(— !)] = V 12 -C'(-l), 
and £(s) is the Riemann zeta function. The latter contribution enters from the integral 

/ t 2 i){t + l)dt = -2 I t\nF(t + l)dt = i(l-ln27r) + 21nA (2.4) 
Jo Jo 2 

This integral may be readily determined from Kummer's Fourier series for In V. Oth- 
erwise, it may be found through the infinite series 

t^{t + l)dt = -1 + £(-l) fc C(*) j( t k ^dt = -| + E YT2 ak) - (2 ' 5) 

k ^ = 2 k — 2 

For reference, we have from (2.1) 
1 1 °° 

TTT, V = « K n + + 2 )( n + 3 K+5 - 3(n + l) 2 (n + 2)p„ +4 

In 1 — z) o ^— ' 

v 7 n=l 

+ (n + l)(3n 2 + 3n + l)p n+3 - n 3 p n+2 ] ^-^q-^+^-^+F' N < L (2 ' 6) 
Then also using (1.3), (1.4), and (2.1), we find 

1 l^ 4 



+ - 



ln(l — z) z / 

= -E[(n+l)(n+2)(n+3K + 5-3(n+l) 2 (n+2)p n+4 +(n+l)(3n 2 +3n+lK +3 -n 3 p n+2 ]^ 

ra=l 

~720 + 2 + + 2)P " +4 " ^ + 1)(2n + 1K+3 + ^"+ 2 ]^ 1 

71=1 

OO OO 

+6^[(« + 2K +4 - (n + lK+3]^"- 1 + 4^p n+ /" 1 . (2.7) 

n=l n=l 

By using (1.2) we obtain 

f 1 ( 1 1\ 4 , f 1 f 1 53 57 2 143 , t 4 

/ —, r + - ) dz = / + — t 1 2 + 1 3 + — 

J \\n(l-z) zj J V 720 18 8 36 24 



+ 



1 8 „ 2 10 o 

1 + Qt 2 1 3 

6 3 3 



[ip(t + l) +7] dt 



49 5 
"72 + 2lRA + 2^ C(3) - 



(2.8) 



From (2.6) we have 



1 1 00 

ln5(1 _ , =^J2^ n+1 ^ n + 2 ^ n + 3 )(™ + 4 K+6 - 2(n + l)(n + 2)(n + 3)(2n + 3)p n+5 

^ ' n=l 

+(n + l)(n + 2)(6n 2 + 12n + 7)p n+4 - (n + l)(2ra + l)(2n 2 + 2n + l)p n+3 + n 4 p n+2 ] z n 

25 



1 5 
+ 



leading to 



24z 8z 2 12z 3 2z A z 



1 

+ - 



51 111 



(2.9) 



ln(l — 2;) 2; 



n=l 



+(n + l)(n + 2)(6n 2 + 12n + 7)p n+4 - (n + l)(2n + l)(2n 2 + 2n + l)p n+3 + raV+2] -2™ 

5 00 

+ -^[(n+l)(n+2)(n+3K +5 -3(n+l) 2 (n+2)^^ 



n=l 



+5 Y^[{n + 2)(n + 3)p n+5 - (ra + 2)(2n + 3)p n+4 + (n + ljV+s]*"" 1 



n=l 



+10 j}(n + 3K +5 - {n + 2) Pn+ ,]z n ~ 1 + b^p n+5 z n -\ 



(2.10) 



n=l 



n=l 



By performing the integration and using (1.2) we determine 



1 1 . • 

+ -) dz 



ln(l — z) z 
Here we have used 



/ 

Jo 



t 4 ij(t)dt 



49 
180 



4367 5. . 15 . /0 , 35 .„ „, 
'8640 + 3 lnA+ 8^ C(3) -T C( - 3) - 



-4C'(-l) + C / (0) + 6C'(-2)-4C'(-3) 



(2.11) 



11 +41nA-iln(27r)-A C(3 )-4C / (-3). (2.12) 



180 2 v ' 2tt 2 



Based upon the evaluation of J Q t k ip(t)dt = —k f Q t k 1 hi T(t)dt, k > —1, we may 
anticipate that J ( ln ^_ z ^ + \ J dz evaluates in terms of 

Q + QlnA + Q^ + Q^ + ... + ^ C(A; " 1) 



7T 2 7T 4 7T fe 2 



and 



+QC / (-3) + QC'(-5) + ... + C'(l-A;), k even, (2.13a) 



+ QlnA + Q^ + Q^ + ... + QM 

7T Z 7T 4 7T K 1 



+QC'(-3)+QC'(-5) + ... + C'(2-A;), k odd. (2.136) 

We provide some discussion of intermediate calculations in the above steps. First, 
we consider the integrals t k ip(t)dt. These may be treated by multiple integrations 
by parts using the fact that —ip(a) is the zeroth Stieltjes constant for the Hurwitz 
zeta function ((s,a). That is, we have the limit representation 

^(f) = -lim(c(z,0--V). (2-14) 
z-»i \ z — 1 / 



Thus, 

' di, (2.15) 



^ t k rl>(t)dt = -limjf * 



A com)- 1 



2 - 1_ 

with the interchange justified by uniform convergence of the integral. We have the 
properties d a ((s,a) = —s((s + l,a) and ((s,t)dt = for Re s < 1. By iteration 
we may then obtain the integrals 

/ t k ((z,t)dt = — / t k d t ((z - l,t)dt 

Jo z - 1 Jo 
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k fl t k -\(z-l,t)dt- ^ Z V (2.16) 



z — 1 Jq ' z — 1 

Secondly, the insertion of the integral representation of (1.2) for p n+ ± into sums 
such as (2.10) gives certain hypergeometric summations. Again, p F q denotes the 
generalized hypergeometric function. As an illustration, consider a contribution from 
a sum such as Y^=i n ^Pn+&- By recalling the property (a) n +i = a(o + l)n, we have 

-E^7=xS = -B» +1 )' H) * M 



J (n + 5)! ^ V 7 (n + 6)! 

M (6 - 0, 



= -t(t-l)(t-2)(*-3)(t-4)(t-5)gi [ j, (n + 6)I 

-t(t - m - 2)( t - 3 )(t - 4)(t - 5)^ f; ^ (6 - t} - 1 



ra=0 



E (7)n n! 



= _ _ t (t - l)(t - 2)(t -3)(t- 4)(t - 5)^ 4^3(2, 2, 2, 6 - t; 1, 1, 7; 1). (2.17) 

Thirdly, (i) there are relations between the p +iF p sums so obtained, and (ii) they 
are often related to the digamma function. By manipulating divided difference forms 
of the tp function, relations such as the following may be obtained: 

3 F 2 (2, 2, 3 - t; 3, 4; 1) + 4 F 3 (2, 2, 2, 3 - t; 1, 3, 4; 1) 
12[3-3* + 7f + t^(t + l)] 12[l-2i + 7* + tV>(t + l)] 12 , 010 , 

H 777 777- — r — 777 77- (.2.I0J 



t(t — l)(t — 2) t(t — l)(t — 2) t(t-l)' 

Likewise, we have 

5 F 4 (2,2,2,2,3-M,1,3,4;1)+ 4 F 3 (2, 2, 2, 3 - t; 1, 3, 4; 1) = ( 2 .1 9 ) 

In this way, summations over the p n constants may be evaluated. 

10 



We also record the following. 
Lemma 3. For Re x > 0, 



~ lux = - J^Pn+iK^- (2-20) 

n=l 



Corollary 2. For Re x > — 1, 



a#(x) - lnx) = - ^Pn +2 , — t • (2.21) 



Proof. For Re x > 0, 



oo / J i 







V»(ac) - In a; = / e"^ ( - - - — — ) dt. (2.22) 



Let t = — ln(l — z), giving 

M x) - iRX= -J\ 1 - z) ->( ] -J-^ + l) dz . (2.23) 

Now use the generating function (1.3), so that 

oo _i 
- lnx = - J^Pn+i / (1 - zf^z^dz 



n=l 



= -$^«+i£(x,n). (2.24) 

n=l 

For Corollary 2, we use x/(x) n = l/(x + l) n -i- □ 
The proof we have given of (2.20) complements that of Proposition 5(a) of [4]. Of 
course this relation may be directly verified with the aid of (1.3): 

v> (» - 12] f'f 



n=l 



11 



-f 

Jo 



(x) - tp(t + x)]dt = ij)(x) - ln[T(x + l)/r(x)] = ^(x) - lux. (2.25) 



We may note that the representation (2.23) may be repeatedly integrated by parts. 
We have for example 



1 f 1 

ip(x) — lax = / (1 — z) 

x Jo 

1 - z)*- 1 

? Jo 



\n A (l-z) z 3 



+ 



1 1 
1 — 

dz 



dz — 



1 



2x 
1 



ln 4 (l - z) z 



6 12 7 1 



z 3 z 2 z 



dz — 



2x 12x 2 
1 1 



2x 12x 2 



x" 



x-l 



24 



+ 



24 60 



ln 5 (l - z) ' z 5 ' ' 



z H 



57 22 1 

z 3 z 2 z 



dz- 



1 



1 



+ 



1 



2x 12x 2 120a; 4 ' 
(2.26) 



In writing these equations, we have used the values of p 2 , ps, Pa, and p 5 for the 
boundary terms. We note that the latter terms give the asymptotic expansion of 
ip(x) — lux as x — > oo. We expand on this point next with regard to Corollary 3. 
The asymptotic expansion 



ip(z)-\nz = -—-^2 



B- 



2n 



1 1 

+ 



2z ^2nz 2n 2z 12z 2 120z 4 252z 6 

n=l 



+ ..., 



(2.27) 



where are the Bernoulli numbers, is well known. In fact, we may readily derive it 
in the following manner. We have 

xjj(x) =#z) -\nx+ = -jT e~ 2te ^cotht- ^ dt. (2.28) 

The asymptotic form of the integral for large x is obtained as t — > 0, in which case 
we may use the expansion 

1 00 2 2k 
coth t = ^ 



R +2fe-l 



+ 2 ^ ^2 
t < 7T . 



(2.29) 
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Then (2.27) follows. Now we may apply Stirling's formula for T to Lemma 2, so that 
for x — > oo, 



^(ar) - lux = - ^p n+l — 

n=l 



nil - n) n(2 - 3n - 2n 2 + 3n 3 ) 
1 H H — - 

2a; 24x 2 



n 2 (— 2 + n + 3n 2 — n 3 — n 4 ) 
-| i i + . . . 

48x 3 



Therefore, by matching asymptotic expansions, we have the following. 



(2.30) 



Corollary 3. (a) p% — A = ^ and ^ = 6p 5 +p 3 — 6p 4 . (b) : y a may be expressed 
as a sum of p n values with rational coefficients. 

Elaborating part (b), we have the following. We let S(n,k) denote the Stirling 
numbers of the second kind. 
Proposition 1. 



^ i = E(- 1 ) fc ( fc - 1 ) ! ^ fc )E^+ 2 

k=l £=0 

n— 1 n 

]>>+ 2 E (-V k (k-l)\S(n,k). 



fc-i 



(2.31) 



£=0 fe=£+l 



Froo/. We let B ( k k) 



-l) k f (x)kdx be the Norlund numbers (e.g., [HI 15]), such 



that Sj 0) = 1, B\ L > = -1/2, and B^> = 5/6, and have 



(2) 



B { : ] + nBt\ l) = (-l) n+1 



nlp n+1 . 



(2.32) 



By iterating, using the initial value B^ = 1, we obtain 



n-l 



BW = (-l)»n! l-J> +2 



(2.33) 



fc=0 
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Substituting into 



and using the sum 



gives the Proposition. 



k=i 



k n 



J2(-l) k (k-l)\S(n,k)=0 



k=l 



(2.34) 



(2.35) 



□ 



The integral representation (2.23) also leads to the following series representation. 
Proposition 2. For Re y > 0, 

1 1 °° 1 

]nF(y) -y\ny + y = -—- -^{y) - J}(ra + 2)p n+3 - n Pn+2 \B(y, n + 1) + - 1h(2tt). 

" n=l 

(2.36) 

Proof. We first integrate (2.23) from x = 1 to y, giving 



lnT(y) -ylny + y- 1 



1 [(i - ^r 1 - 1] 



ln(l — z) \ m (l ~ z ) z 



dz 



- / [(1 - z)y- 1 -1] l)p„ +3 - np n+2 ]z n +p 3 -- + Vp n+1 /- 2 } dz 

J ° U=l ^ n=l J 

„1 f oo oo ^1 

= - / [(l-z)y- 1 -l}{J2[(n + l)Pn + s-n Pn+2 }z n +p 3 -- + J2Pn + 2Z n - 1 \dz : 

J U=l ^ n=l J 

(2.37) 

using (1.3) and (1.4). Performing the integration, we find 



lnr(y)-ylny+y-l = - ^[(n+l)p n+3 -np n+2 ] 



n=l 



S(y,n + 1)- 



n+1 



-P3 



n=l 



5(y,n) - 



(2.38) 
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Next, the integral representation (1.2) is used for the sums absent the Beta function, 

£^ = ^(270-1-7, (2.39) 



Pn+2 _ 1 

n ~ 2 

n=l 

5 n 



£pn+3 = ^. "£^Y^+2=7-l- (2-396) 

n=l ~" ra=l 

Alternatively, these sums may be obtained by integrating and other otherwise ma- 
nipulating the generating function (1.3) and taking z — > 1. Then combining terms of 
(2.38) and (2.39) gives the Proposition. □ 
The coefficients p n +i may be readily related to other quantities, including the 
Bernoulli numbers of the second kind b n (n > 0) [H [TO] . 

b n = I r ( t + 1 ) dt, (2.40) 
7 T(t-n+l) ' K J 

with &o = 1, &i = 1/2, &2 = -1/6, and b 3 = 1/4. We have 

1 /' 1 r(n-t), (-l) n+1 Z" 1 r(t + l) , , , , 

Pn+i = — r / —^-—^-dt = / V ; dt, 2.41 

and hence p n +i = (— l) n_1 &n/^!- 

The following gives a series representation for In A. The method of proof shows 
that a family of such series may be written. 
Proposition 3. 

, , 1 AfL x , \ n 2 (6n+ 1) 1 1 

\nA = -+y < -(n + 2)p n+i + (n + 1) 1 »n+3 + — ? r - 1 + —, t Pn+2 ^ 

4 f^l 2 ^"+ 4 ^v -r y ^ n+3 ^ [2(n + 1) 12(n+l)J^ + J 

(2.42) 
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Proof. Let C( s > a ) be the Hurwitz zeta function. Then for Re s > —(2n — 1), 
n G N , and Re a > 0, there is the integral representation 



B 



((s,a) = a- s + J2(s)k-ijla- k - s + 1 + 



k=0 



\s)Jo \e*-l ) 



By taking n = 2 and a = 1, one may find 

1 



1 * 



1 Z" 00 

In A = - + / . , 

4 ./n le'-l t 2 12 7 t 2 



(2.43) 



(2.44) 



Changing variable with t = — ln(l — z) we have 



1 



\nA = -+ (- + 

4 7 ^ ln(l-z) 2 



_ l + ln(l-*) 



(2.45) 



12 ;in 2 (l-z)' 

Next the generating functions (1.3), (1.4), and (2.1) are employed in the integrand, 
and all terms 0(z~ k ), k = 1,2,3 are nullified, as they should. Performing the inte- 
gration gives 



71=1 

11=1 
1 oc 



n + 1 



n 



£>n+3 - Pn+2 



n 



n + 2)p n+4 - (2n + l)p n +3 H — Pn+2 

n + 1 



71=1 



Pn+3 



n+ 1 



+ P3 



1 



Pn+1 



12 ^ n 

n=l 



(2.46) 



Shifting the index on the last sum and combining terms gives the Proposition. □ 
Noting that In A ~ 0.248754477033784262547253, the summation in (2.29) pro- 
vides an appropriate correction to 1/4. 
Corollary 4. The constant 



C'(2) = C(2)[ 7 + ln(27r)-121nA] 



(2.47) 



16 



may be written in terms of the coefficients p n . 

Proof. The relation of C'(2) to ('(—1) follows from the functional equation of 
the Riemann zeta function. The value £(2) may be found in terms of p n 's via the 
integral representations (2.34) at a = 1. The constant ln(27r) may be written in 
terms of p n 's via (1.7) and (1.9). Finally, the Euler constant 7 = Y^=i "^Ti as 



We may note that more generally we may similarly write series representations 
for the logarithm of the double Gamma function T 2 , since we have for Re a > 







1 

- + 

z 



ln(l - z) 2 



1 1 




0-1 



ln 2 (l-^) 



dz 



(2.48) 



We arrive at 



Proposition 4. For Re a > 0, 





(2.49) 



The details of this calculation are omitted. 
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We conjecture that the following inequalities hold for the coefficients p n . 
Conjecture 1. (i) For n > 2, 

(n + l)p n+3 - np n+2 < 0, (2.50) 

(ii) that the sequence {p n } is strictly convex, i.e., for n > 3, 

Pn < 7>(Pn+l +Pn-l), (2.51) 

(iii) that the sequence {p n } is strictly log-convex, i.e., for n > 3, 

p 2 n < p n „ip n+1 . (2.52) 

The convexity and log-convexity themselves are probably not difficult to show, and 
(iii) implies (ii). Certainly these inequalities hold for n sufficiently large. Using 
the known asymptotic form p n ~ l/[n(lnn + 7) 2 ], the leading asymptotic form of the 
difference of the left and right sides of the Conjecture is given by: for (i), — 2/[n(lnn + 
7) 3 ], for (ii), — In 2 n/[n 3 (Inn + 7) 4 ], and for (iii), — In 2 n/[n 4 (Inn + 7) 6 ]. 

Expressions for the Stieltjes constants and ((m,a) values 

The Hurwitz zeta function, defined by ((s, a) = Yln^oi 72 + a )~ s f° r Re s > 1 and 
Re a > extends to a meromorphic function in the entire complex s-plane. In the 
Laurent expansion 

C(s, a) = - + [ -^^(a)(s - If, (3.1) 

s — 1 L — ' nl 

n=0 

7„(a) are the Stieltjes constants (e.g., OH]), and by convention one takes 7^ = 7fc(l). 
One has 70(a) = —^(a) and 70 = 7. The Stieltjes constants may be expressed via 
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the limit formula 



7fc(a) = lim 

Af-5>oo 



j^\n k (n + a) ln k+1 (N + a) 



n=0 



71 + a 



We dispense with further preliminaries concerning the 7 n (a)'s. 

We may write series and integral representations for these constants based upon 
the p n coefficients. This development is illustrated in the next result. 
Proposition 5. (a) 



- ln 2 a- 7l (a) = 7 [ln a - ^(a)] + $>n+i / « a_1 ln(- lnu)(l - m)™" 1 

~ n=l ^° 

= 7 [In a - V'(a)] + / u a ~ l ln(- ln-u 



1 1 

1 — u ln-u 



and 



-7i-7 



/ e-*lnt(l-e-*)^ 1 ^ 

n=l J ° 



POO 


- 1 e-*" 






Jo 


1 -e* + t _ 



lnt rft, 



(b) 



and (c) 



71=1 



2 7_! 
2 



-73 = 7 4 + y7i + — ( 71_2 + 6 7i) + 3 772 + 27((3) 



/>CO 00 

+ / (l^^e-^^+^l-e-*)™- 1 ^. 
Proof. For Re s > and Re a > we have the integral representation 



C(s,a) - 



(s - 1) a- 1 r( s ) 7 



1-e" 



— ) c/.c 



(3.2a) 



(3.26) 



72 = -7(7 2 + C(2) + 2 7l )+ / (ln^e-'^p^a-e-T- 1 ^, (3.3) 



(3.4) 
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, oo „i 

= ^Vp n+ i ^(-ln^-^l-ur-Mu. (3.5) 
1 W „=i 

It is readily found that the limit of (3.5) as s — > 1 agrees with Lemma 2, and 
this result is used in all parts of the Proposition. For parts (a)-(c) we take succes- 
sive derivatives of (3.5) with respect to s and evaluate at s = 1. For the integral 
representations in (a) we use the integral representation for p n+ i in (1.2). □ 

The following further emphasizes the connection of the p n+ i coefficients with an- 
alytic number theory. We let ip' be the trigamma function, ip^ the polygamma 
functions, and H n = Ylk=i \ = "^( n + 1) + 7 be the nth harmonic number. 
Corollary 5. 

oo 

C(2) = l + ]T^±itf n , (3.6) 

n=l 

and 

C(3) = \ + ^7 + \ E P -^[H 2 n - rl/(n + 1)]. (3.7) 
2 2 2 z — ' ra 

n=l 

Proof. This follows from the special case of (3.5), 

1 1 00 Z" 1 

C(*)"(^= j^£Pn+i j( (-Inur^l-u)- 1 ^. (3.8) 

The integrals are given by 

/ {-lnu) a - 1 {l-u) n - 1 du= {-ly^d*- 1 B{x,n)\ x=1 . (3.9) 
Jo 

In writing (3.7) we have used the previously given summation expession for 7. □ 



1 

- + 



z \n(l-z) 



dz 
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Similarly further values of C(n) may be written. We may note the appearance of 
generalized harmonic numbers = Ylk=i W-> ^« = Hrp . In particular, as regards 
(3.7), HL 2) = -[f{n + 1) - V'(l)] = -W(n + 1) - C(2)], so that 



1 , 1 Pn+1 

<,w = 

Furthermore, 



C(3) = ^£^[tfn + tff]- (3-10) 

n=l 



n=l 



As in Corollary 5, the summand is positive and the convergence is from below. 
The generalized harmonic numbers are given by 



HP = [^-D(n + 1) -^(l)] . (3.12) 

We may systematize Corollary 5 and related results in the following manner. We 
introduce the (exponential) complete Bell polynomials Y n = Y n (xi, x 2 , ■ ■ ■ , x n ) ap- 
pearing in the expansion 



CO 

exp 



y^x m —)=l + y^Y n (x 1 ,x 2 ,...,x n ) — , Y = l. (3.13) 

m=l / n=l 



Then 

Proposition 6. For m G N, m > 2, (a) 

^ ' ■ n=l 

(3.14) 

and (b) for Re a > 0, 

C(m,a) = ^^p n +iB(a,n)Y m _ 1 [^(a)-^(a+n),'ip'(a)-'ip(a+n), 

[m — 1) a m 1 (m — \)\ 
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^ m ~ l \a) - V (m_1) (o + ri)]. (3.15) 

Proof. We apply (e.g., [6]) 
Lemma 4. For differentiable functions / and g such that f'(x) = f(x)g(x), assuming 
all higher order derivatives exist, we have 

f {x) = f {x)Yj [g(x),g'(x),...,g^(x)] . (3.16) 

We put f(x) = B(x, n) and g(x) = In -B(x, n) = ip(x) — ip(x + n). Then by (3.12) 
g ( r )(l) = (-l) r (r - l)\H { n\ Part (a) then follows from (3.8) and (3.9) with s = m 
and -6(1, n) = 1/n. 

Similarly for part (b) we use (3.5), in which the integral on the right side is given 
by {-ly^d^Bfan). □ 

Alternatively, the integral on the right side of (3.5) may be treated with a gener- 
ating function for the Stirling numbers of the first kind [1] (p. 824) so that 

/ (l-u) a - 1 [-\n(l-u)] m - 1 u n ~ 1 du= (-l) m -\m-l)\ Y] ±-L s (j, m -l)B(n+j,a) 

J ° j=m-l J - 

(3.17) 

Additional sums 

We collect the following summations with the p n +i coefficients. 
Proposition 7. 



l'n + \ I 

n=l 

in particular 



V = - - 1 1 B(a, x + l)dx, (3.18) 

z — i n + a a 

n=l 







yPn ±L = 1 _ ln2 (319) 
^ n + 1 v 7 



„ , ■- + 
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n=l 

oo 



00 1 2 1 /-l 



™ 3 12 QJo 



n=l 

and for j G N 



V = z / x j+2 F j+1 (l, 1, . . . , 1, 1 - x; 2, 2, . . . , 2; z)dz. (3.22) 

n=l ™ ^ 

Proof. These may be obtained with the aid of the integral representation of (1.2), 
and we omit further details. □ 

As concerns the integral on the right side of (3.20), we have the following result. 
Corollary 6. 



1 ^2 r\ 

1 ■ \ 1 1 1 

dx. (3.23) 



i) 2 (x + l)dx = 271 - — + 



2^ + \-2^(x)+7f/(x) 

X X A 



Proof. We use Proposition 3(b) of [5], 

7i = %■ + / (tI>{x) + - il>'{x)] ) dx, (3.24) 



together with i/j 2 (x + 1) = V> 2 (z) + 2^ + 4j. □ 
Based upon various representations of the Beta function, we may collect the fol- 
lowing representations for the integral occurring in (3.5) for the Riemann zeta function 
case, Qj = 1. It should be clear which of these are restricted to s > 1 an integer. 
Proposition 8. 

i n ( s )= [\-\n U y-\i-ur- i du= (-lr-^r^MU 

Jo 



2.3 



T(s) >A (~l) k - 1 fn 

2-^ hs-i 



n * — ' k & 

k=l 



n *■ — ' 7i • • • 7 s -i 



r(s) / (1 — Xix 2 • • • x s ) n l dx\---dx s 

'[0,l] s 



^r>) 

2ni n 



1 n 

in 



dz 



\z\=r<l Z S fJl (1 - Zfj) 



1/2+ioo 



;-i) w r( a ) 

2tt2 n A/2-ioo 2/ s (2/ - 1) • • " (y - n) 



dy. 



(3.25) 



Furthermore, for s > 2, 



n 



s - 1 
n — 1 



+ - l), 
n 



and 



n 



(3.26a) 



(3.266) 



Upper bound for p n +\ and other results, including 
polylogarithmic representation of the Stieltjes constants 

We let Si(z) = f* ^dt be the sine integral. 
Proposition 9. For n > 1, 



Pn+l <-\ + -Si(vr) + (-1)- 1 

L 7T 



1 Si[(n - 1)tt] 
2 



7T 



(4.1) 



Proof. Based upon contour integration, Knessl has shown the following integral 
representation [3], 

1 du 



Pn+l 



o (l + w) n (ln 2 w + 7r 2 ) 
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n > I. 



(4.2) 



Then 



f 1 1 du f 00 1 du 

Pn+1 ~J (l + u) n (ln 2 M + 7r 2 ) + A (l + M)"(ln 2 « + 7r 2 



< 



du 



+ 



(lll 2 M + 7T 2 ) J 1 u n {ln 2 u + 7i 2 ) 

du 



"00 , J J 



f 

V^ 2 « n / (ln 2 M + vr 2 ) 
= Jo v 2 + ^ 2 



-dv 



(4.3) 



= ~ + ^Si(7T) + (-l) 
Z 71 



n-1 



1 Si[(n - l)?r] 
2 



7T 



The integral of (4.3) could be evaluated by means of contour integration, but we 
supply another means. We put 

»oo e _(fc_l)„ 



and form 



/(*) 



/"(ife) + 7l 2 I(k) 



V 2 + 7T- 



e-^dv = 



1 



fc - 1 



(4.4) 



(4.5) 



The homogeneous solutions of this differential equation are of course cos kn and sin kn 
and their constant Wronskian is W = n. Per variation of parameters, a particular 
solution then takes the form 

T ,,s cos7r/c /" sin7r/c sin7r/c /" cos7r/c „ 
ip(A;) = / -dk H / -dk. 



Tt 



k-l 



Tt 



k-l 



Then solving (4.5) subject to 1(1) = 1/2 gives 



I(k) = — {— sin(A;7r)Ci[(A; — l)n] + cos(A;7r)Si[(A; — l)n]} — - cos kn + c 2 sin &7T, (4.6) 

7T 2 
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where Ci(z) = — f°° ^^dt is the cosine integral. Then imposing I(oo) = gives the 

constant C2 = 0. □ 

Sif 

Asymptotically, as n — > 00 on the right side of (4.1), this upper bound is — — 

We may extend Lemma 3 to the following. 
Lemma 5. 



x x) n 2x J (t 2 + x 2 )(e 27rt - 1) 

-^_2 r °° Vdv 



2x J (1 +t» 2 )(e 2 ™ - 1) 
I r 00 lnu 



2x A u(l + ln 2 u)(u 2 ™- 1) 
_1 Z" 00 2 F 1 (l,l;x + l;v)dv 
xj W [ln 2 (i-l)+7r 2 ]- l4 '° 

Proof. The first line follows from a known integral representation of ip(x) — In a; 
[7] (p. 943). The final equality follows from the representation (4.2) for p n+ %. □ 

For the next result we introduce the polylogarithm function Li a (z) = Y^kLi f? 
that may be analytically continued to the whole complex plane. This function has 
the integral representation for Re s > 

u ^=w)L —• (48) 

and we note that IAi(z) — — ln(l — z). The branch cut for IA s (z) in the complex z 
plane may be taken from 1 to 00. This function is also given when s is a positive 
integer by Li m (z) = z m+ iF m (l, 1, . . . , 1; 2, . . . , 2; z), and this provides another way of 
seeing the branch point at z — 1. 
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Proposition 10. 



and 



vi - r = I " in ( i • - ) • d s \ s=1 u s 



72 + 7[7 2 + C(2) + 2 7l ] 



i : r^2))\n[l + -) +2 1 d s \ s=1 U s (-- 

Hi J \ Uj 



Proof. By (3.2b) and the representation (4.2), 



du 



(ln 2 « + 7r 2 )' 



(4.9) 



di\_. Li s — 
1 v u 



du 



(ln 2 M + 7r 2 ) 
(4.10) 



-li-l 2 = Y ; s 7 2 7/ e-*lni(l-e-*r -1 ^ 



du 



In t (it. 



/o Jo (! + e * M ) 

The interchanges are justified by the absolute convergence of the integrals. Then (4.9) 
follows from logarithmic differentiation of the integral of (4.8). Similarly for (4.10), 



(3.3) is used along with (4.2). 

Several versions of (4.10) may be written by employing (4.9). 

Proposition 10 is clarified and generalized with the following. 
Proposition 11. We have 

1 f°° U s (-v)dv 



□ 



Consequently, for k > 0, 

Ik = (-1) 



S — 1 Jo f 2 (ln 2 V + 7T 2 



(4.11) 



k-ifd 



ds 



s=l 
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Li s (— v)dv 
o f 2 (ln 2 v + 7T 2 ) 



(4.12) 



Proof. We first combine the representation (3.5) with a = 1 with the second line 
of Proposition 8 resulting from binomial expansion: 

Pn+i A fn\ (-l)^ 1 



n=l fc=l v 7 



Pn +2 yfn + iN (-i) fc >p ^MHf 

^ (n + 1) ^ U + 17 (A; + l) 5 " 1 2^P n + 2 2^\k)(k+ IV 

n=0 K ' k=0 V 7 V ' n=0 fc=0 V 7 V ' 



OO OO 



k=0 n=k \ / \ I 

We now use the representation (4.2) for p n+ 2, so that 



w.x _ j_ = v ( - 1)fc 1 du = _ r Lis ^ 

UJ s-1 u k + l (k + I)* (In 2 U + TT 2 ) Jo 



k -< o u k+l (k + l) s (ln 2 M + 7r 2 ) Jo ln 2 u + n 2 
Li s (— v)dv 



du 



/o f (hi V + 7T Z ) 

(4.12) then immediately follows. □ 
Example. When k = in (4.11) and s — y 1, 



» = ^ 1+ E> . < 4 ' 15 ' 

f Z (m t> + 7T Z > 



and this case is equivalent to (2.87) in [3]. 

Let $ denote the Lerch zeta function, s, a) = Yl^Lo ( n +a) 3 ' anc ^ analytically 
continued. This series holds for s G C when |z| < 1 and for Re s > 1 when \z\ = 1. 
The function $ has an integral representation 

1 f°° t s ^e- at r 1 ^-i e -(«-i)* 
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for Re s > when \z\ < 1, 2^1, and for Re s > 1 when z — 1. 
Proposition 12. 



. / x 11 Z" 00 $(-u,s,a) , 

C(s, a) - r = / n 2 ' > • ( 4 -!6) 

(s — lja s 1 y v(m v + rr 2 ) 

Proof sketch. Now 



Then (3.5) is used, following steps similar to the proof of Proposition 11, so that 



v ' fc=0 n=fc \ / V / 

Performing the sum on and changing variable in the integral gives (4.16). □ 
Corollary 7. 

In a — ^(a) = In a + 70(a) = / — - — 5— ^ — , ——dv. (4.19) 

Jo v(\n v + 7T 2 ) 

Appendix 

Here we show that for n G N + 

(a)j _ (a) n ^fn-l\ (-l) k (a + k), 



( a + n ). ( n _i)|^^ k J( a + k)(a + k + l) J 
Proof. We have the known decomposition 

1 1 » ( N \ (-ir 

x(x + l)---(x + N) N\f^\k J (x + k)' K ' J 

Then 



(a)j a(a + l)(a + 2)---(a + n-l) 



(a + n)j (a + j)(a + j + 1) • • • (a + j + n - 1) 
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(q)n y^(n-i\ (~i) fc U3) 



By noting that 

(a + k)j (a + fc) 



(A4) 



{a + k + l)j (a + k + jY 
the result follows. □ 
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